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Abstract 

It is proved that, in two dimensions, the Calderon inverse conduc- 
tivity problem in Lipschitz domains is stable in the L p sense when the 
conductivities are uniformly bounded in any fractional Sobolev space 
W a > p a > 0, 1< p < oo. 
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1 Introduction 

Calderon inverse problem consists in the determination of an isotropic L°° 
conductivity coefficient 7 on 0, from boundary measurements. These mea- 
surements are given by the Dirichlet to Neumann map A 7 , denned for a 
j. ■ function / on dQ as the Neumann value 



o 
00 
o 



A 7(/) = 7 — U, 



d_ 
du 



where u is the solution of the Dirichlet boundary value problem 

|v-( 7 Vn)=0 
1 u\on = f 



and ^j denotes the outer normal derivative. For general domain and conduc- 
tivities where the pointwise definition 7^it has no meaning, the Dirichlet 
to Neumann map 

A 7 : H^ 2 (dn) -» H-V 2 (dn) (1.2) 

can be defined by 

(A 7 (/), V0 )= / 7Vn-Vv9 (1.3) 

where ip € IF 1,2 (f2) is a function such that ip\gn = ipo in the sense of traces. 



Since the foundational work of Calderon, the research of the question 
has been very intense but it is not until 2006 when, by means of quasicon- 
formal mappings, K. Astala and L. Paivarinta in [llj were able to establish 
the injectivity of the map 

7 ^A 7 

for an arbitrary L°° function bounded away from zero. Previous planar 
results were obtained in [39] and |3S]. In higher dimensions, the known 
results on uniqueness require some extra a priori regularity on 7 (basically 
some control on | derivatives of 7, see [17] , [15] . [IT] and [2"T].) 

A relevant question (specially in applications) is the stability of the 
inverse problem, that is, the continuity of the inverse map 

A 7 ^ 7 . 

For dimension n > 2, the known results are due to Alessandrini [1], [5]. 
There the author proved stability under the extra assumption 7 € W 2 ' 00 . 
In the planar case, n = 2, the situation is different. Liu proved stability 
for conductivities in W 2 ' p with p > 1 in [36] . In [Jl], stability was obtained 
when 7 G C 1+a with a > 0. Recently, Barcelo, Faraco and Ruiz [15] obtained 
stability under the weaker assumption 7€C a ,0<a<l. Precisely, they 
prove that for any two conductivities 71,72 on a Lipschitz domain f2, with a 
priori bounds -^ < 7« < AT, A' > 1 and ||7j||c q < Ao> the following estimate 
holds: 

II71 -12\\l°°{0.) < ^(||A 7l - A 72 || i/ i/2( 9n) _ >jff -i/2( 9n) ) 

with T/(t) = C log(|) _a . Here C, a > depend only on if, a and Ao. 

An example, due to Alessandrini [I] , shows that in absence of continuity, L°° 
estimates do not hold. Namely, if we denote by B ro = {x G R 2 , \x\ < r^} 
the ball centered at the origin with radius vq, take Q = B\ the unit 
ball in R 2 , 71 = 1 and 72 = ■■ 1 + XB Tn , then H71 - ^h^m '- : 1, but 
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|A 7l — A 72 || rr i rr _i < 2ro — > as ro — > 0. 



A closer look to the previous example shows that lim ro ^o ||7i — T2||z, 2 (si) = 0- 
Therefore one could conjecture that, in absence of continuity, average sta- 
bility (in the L 2 sense) might hold. However, it is well known that some 
control on the oscillation of 7 is needed to obtain stability. Namely, let 7 
be defined in the unit square and extended periodically. Then the sequence 
{ i y{jx)}^2- 1 G-converges to a matrix 70 (see for example [50] for the notion 
of G-convergence) . On one hand, j(jx) has not any convergent subsequence 
in L 2 . On the other hand, G-convergence implies the convergence of the 
fluxes [50} Proposition 9]. That is, if Uj,uq solve the corresponding Dirichlet 



problems for a fixed function / £ H? (dfl), 

(1.4) 



v • ( 7j -VttjO 



u 



,jl9fi = / 



then, the fluxes satisfy that jjUj — *■ 7VU. Thus, by (|1.3p 

limj 1 j 2 _ >00 (7V 7i — A 7j ),v?o) for each yjQ- Notice that jj can be cho- 
sen even being C°°, so the problem here is not so much a matter of 
regularity but rather a control on the oscillation. 

In this paper we prove that L 2 stability holds if we prescribe a bound 
of 7 in any fractional Sobolev space W a ' 2 . By the relation with Besov 
spaces this could be interpreted as controlling the average oscillation of the 
function. Thus average control on the oscillation of the coefficients yields 
average stability of the inverse problem. 

Theorem 1.1. Let be a Lipschitz domain in the plane. Let 7 = 71,72 be 
two planar conductivities in Q satisfying 

• (I) Ellipticity: ± < 7(3;) < K. 

• (II) Sobolev regularity: 7$ € W a ' p (Q) with a > 0, 1 < p < oo ; and 

\hi\\w a 'P(n) < r . 

Let a = min{a, |}. Then there exists two constants c(K,p), 

C(K,a,p,To) > 0, such that: 

ll7i-72ll^ffi)< |log(p)|cfi2 (1-5) 

where p = ||A 7l - A 72 ||^i/2 (9Q) ^ H -i/2 (0n) . 

The theorem is specially interesting for a — ► 0. Then we are close to get 
stability for conductivities in L°°. An estimate for the behavior of C in 
terms of To is obtained in the particular case p = 2 (see Corollary 16. 3D . 
and analogous results for p / 2 can be deduced by interpolation. Also 
interpolating one can obtain L p stability estimates, whose behavior in a 
will be quadratic as well. 

Concerning the logarithmic modulus of continuity, the arguments of 
Mandache [38| can be adapted to the I? setting. Namely we can consider 
the same set of conductivities with the obvious replacement of the C m 
function by a normalized W a ' 2 function. The argument shows the existence 
of two conductivities such that H71 — 72||l°°(b) < e > Il7i|lvy a -p(n) < To; but 

Il7i - 72||l2 (d) > aa+ii ■ ( L6 ) 

C|log(p)| 2q 



Here C is a constant depending on all the parameters. Notice that the power 
is better than in the L°° setting but still the modulus of continuity is far 
from being satisfactory. 

In our way to prove Theorem 1 1.1 1 we have dealt with several questions related 
to quasiconformal mappings of independent interest. More precisely, we have 
needed to understand how quasiconformal mappings interact with fractional 
Sobolev spaces. In particular we analyze the regularity of Beltrami equations 
with Sobolev bounds on the coefficients which has been a recent topic of 
interest in the theory. See [Ml ES] where the case /j, E W 1,p is investigated 
in relation with the size of removable sets. We prove the following regularity 
result. 

Theorem 1.2. Let a € (0, 1), and suppose that /i, v £ W a ' 2 (C) are Beltrami 
coefficients, compactly supported in D, such that 

IIm(*)I + K*)II< f^j- 

at almost every z € B. Let (p : C — > C be the only homeomorphism satisfying 

d(j) = n dcj) + v dcf) 

and 4>(z) — z = 0(1/ z) as \z\ — » oo. Then, 4>(z) — z belongs to W l+Sa,2 (C) 
for every 9 G (0, i), and 

\\D 1+ea (cf) - z)|| L 2 (c) < C K (j|H!iW,2( C ) + \W\\w^(c)) 
for some constant Ck depending only on K . 

Many corolaries can be obtained form this theorem by interpolation, as for 
example what do you obtain if /x is a function of bounded variation. We 
have contented ourselves with the L? setting but similar results hold in LP. 
As a consequence of this theorem, we obtain the corresponding regularity of 
the complex geometric optics solutions. 

The other crucial ingredient in our proof is the regularity of \i o ip where ij) is 
a normalized quasiconformal mapping. It is well known that quasiconformal 
mappings preserve BMO and W 1 ' 2 but it is not clear what happens with 
the intermediate spaces. We prove the following stament, 

H g W a ' 2 => fio^pe W?' 2 , for every < -^ (1.7) 

K 

which suffices for our purposes. The proof relies on the fact that Jacobians 
of quasiconformal mappings are Muckenhoupt weights [10] 

The Lipschitz regularity of the domain Q is used to reduce the prob- 
lem to the unit disk D. This reduction relies on two facts. First, any 



Lispchitz domain £1 is an extension domain for fractional Sobolev spaces. 
Secondly, the characteristic function xn belongs to W a,2 (C) for any a < ^- 
Indeed, this is responsible also of the constraint a < \ at Theorem ll.il In 
fact, a stability result holds as well if Q is any simply connected extension 
domain. To see this, recall that planar simply connected extension domains 
Q are quasidisks ([28]), that is, f2 = </>(B) where (j) : C — > C is quasicon- 
formal. Therefore, for instance by our results in Section [H \£i = Xo ° 4 > ^ 1 
belongs to some space W a ' 2 , and then use Theorem ll.il 

The rest of the paper is organized as follows. In Section [2] we recall 
previous facts from [11| [T5] which will be needed in the present paper, 
and describe the strategy of our proof. In Section [3] we reduce the 
problem to conductivities 7 such that 7 — 1 G Wq' (ID). In Section 0] we 
study the interaction between quasiconformal mappings and fractional 
Sobolev spaces. Finally in Section [5] we prove the subexponetial growth of 
the complex geometric optic solutions and in Section[6]we prove the theorem. 

In closing we remark several issues raised by our work. The first one 
is to improve the logarithmic character of the stability. It was proved by 
Alesssandrini and Vesella that often a logarithmic estimate yields Lipschitz 
stability for some finite dimensional spaces of conductivities. However, to 
achieve the desired estimates in our setting seems to require a more subtle 
understanding of the Beltrami equation and we leave it for the future. It 
will also be desirable to obtain LP estimates in terms of W a ' p with constants 
independent of p, so that the C a situation in [15] could be understood as a 
limit of this paper. This seems to require an L 2 version of the boundary 
recovery results of Alessandrini [5] and Brown (see [19]). Finally, from the 
quasiconformal point of view, there seems to be room for improvement 
in our estimates specially concerning the composition which is far from 
being optimal when a /" 1, since W 1 ' 2 is invariant under composition with 
quasiconformal maps. This will also be the issue for further investigations. 

Notation 

Complex and real derivatives are denoted by 

dz = d = — = -(— + i— 

&z 2 \dx dy 



d z = d 



d 1 / d . d 



1- 



dz 2 \dx dy 



where z = x + iy. For a mapping (p : Q — > C, its Jacobian determinant 
is denoted by J(z,(p) = \d z (j)(z)\ 2 — \dz4>(z)\ 2 . For k € C we will use the 
unimodular function eu{z) = e tkz + lkz . Notice that then we can define the 



Fourier transform by 



f(k)= / e. k (z)f(z)dA(z). 

Jc 



The spaces 17(0), W^(Q) and W^ip) are denned as usually. Then, 
following Adams pQ, one introduces W a ' p (Q) as the complex interpolation 
space 

w a > p (n) = [L p (n),w 1 > p (n)] a , 

and similarly for the homogeneous case W a,p (Q) = [L p (ft), W 1,p (Q)] a . The 
Holder space C Q (Q) over a domain ft is 

C*(0) = (/ : [|/[|l- + sup l/( , X) "{y < 00} . 
[ x, y eU \x-y\ a J 

For simplicity, H l (ft) = W 1 ' 2 ^) and H% = Wq' 2 (Q). By H^(dtt) we de- 
note the quotient space H^-^ty/H^fl). Given a Banach space X we denote 
the operator norm of T: X — » X by ||T||x- We remark that C or a denote 
constants which may change at each occurrence. We will indicate the depen- 
dence of the constants on parameters K, T, etc, by writing C = C(K, T, ...). 
Finally, for two conductivities 71 and 72, we write 
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2 Scheme of the proof 

We will follow the strategy of |T5] . This work focusses on the approach based 
on the Beltrami equation iniciated in [11] . The starting point is the answer 
to Calderon conjecture in the plane obatained by Astala and Paivarinta. 

Theorem 2.1 (Astala-Paivarinta). Let $7 C M. 2 be a bounded simply con- 
nected domain, and let 7^ £ L°°(Q), i = 1,2. Suppose that there exist a 
constant K > 1 such that -^ < 7^ < K . If 

A71 = ^72 



then 71 = 72 • 

In other words, the mapping 7 1— > A 7 is injective. We recall the basic 
elements from [IT] needed in the sequel, also the strategies for uniqueness 
and stability, and what we will need in the current paper. 

Equivalence between Beltrami and conductivity equation: Let D 

be the unit disc. If a function u is 7-harmonic in D, then there exists another 
function v, called its 7-harmonic conjugate (and actually 7 _1 -harmonic in 
O), unique modulo constants, such that / = u + iv satisfies the M-linear 
Beltrami type equation 

df = fxdf (2.1) 

with 

1-7 
a = - G R. (2.2) 

1 + 7 
Then if K > 1 is the ellipticity constant of 7 we denote by 

K-l 



K + l 



It is an algebraic fact to show that ||/z||oo 5; K an d thus the Beltrami equation 
is elliptic when so is the conductivity equation and viceversa. Moreover, for 
x € {^,K), the function F(x) = ^f satisfies -^ < \F'(x)\ < ^ R . Thus, 
it also follows that 

^7 ll7lliy a 'P(n) < ||/i|lw».*'(ft) ^ C ll7llw«'P(n)) 

where the constant C only depends on K (see Lemma 13. ip . Therefore, 

bounds in terms of fi and 7 are equivalent. 

We can argue as well in the reverse direction. If / € W,' (B) satisfies (|2.ip 

for real fj, with ||/i||oo < K > then we can write / = u + iv where u and v 

satisfy 

div (7 Vit) = and div (7 _1 S/v ) = 0. 

Thus, it is equivalent to determine either 7 or fi, and throughout the paper 
we will work with either of them indistinctly 

As for holomorphic functions, u and v are related by the correspond- 
ing Hilbert transform 



H,,: H*(dB) -» H2(dl 



defined as 



Ti-fiiulan) = v\go 

for real functions, and R- linearly extended to C- valued functions by setting 
H/x^iu) = iTi-^u). Since Bt'H^ = A 7 it follows [11, Proposition 2.7] that 



Hfj,, Ti-fj, and A 7 -i are uniquely determined by A 7 . Accordingly in [15| 
Proposition 2.2] it is shown that 

ll'Ty n__/ II < II A A II 

Il'l/Xl ''•/I2II r^j ll iV 7l '' i 72ll 

with respect to the corresponding operator norms. In other words, the 
mapping A 7 1— > TL^ is Lipschitz continuous independently of the regularity 

of 7. 

Existence of complex geometric optics solutions, scattering trans- 
form and dk equations: The theory of quasiconformal mappings and 
Beltrami operators allow to combine in an efficient way ideas from complex 
analysis, singular integral operators and degree arguments to prove the ex- 
istence of complex geometric optics solutions with no assumptions on the 
coefficients. 

Theorem 2.2. Let k € (0,1), and let /x be a real Beltrami coefficient, 
compactly supported in B, satisfying ||/i||oo < K - For every k € C and 
p € (2, 1 H — ) the equation 

df = fxdf 

admits a unique solution f £ W,'^(C) of the form 



such that M[t(z, k) — 1 



and ffj,(z,0) = 1. 

In this context, the proper definition of scattering transform of \i (or of 7) 
is 

*>.(*) = £ J §~ z (^ Z (U® - f^M) dA(z). (2.3) 

The complex geometric optics solutions {u 7 , u 7 } to the divergence type equa- 
tion (jl.ip are then obtained from the corresponding ones from the Beltrami 
equation by 

Uj =Re(/ M ) +zIm(/_ jU ) 

Uj = Im(/ M ) + zRe(/_ M ), 

and they uniquely determine the pair {/ M ,/- M } (and viceversa) in a stable 
way. We consider u 7 as a function of (z, k). In the z plane, u 7 satisfies the 
complex 7-harmonic equation, 

div(7Vu 7 ) = 0. 



J C VV lOi 


; V^J U J 




jc/ //£ 


/(*) = 


-- e ikz M 


«(*> 


fc) 


0(l/z) 


as \z\ — s 


' 00 


. Moreover, 


Re( 




>( 


) 



As a function of fc, n 7 is a solution to the following d- type equation 

-£L(z,k) = -iT f ,(k)u(z,k). (2.4) 

ok 

Let us emphasize that r M (/c) is independent of z. 

Strategy for uniqueness: Let 71,72 be two conductivities. In |llj . the 
strategy for uniqueness is divided in the following steps: 

(i) Reduction to ID. 

(ii) If A 7l = A 72 , then r Ml = r m . 

(Hi) Step (ii) and (|2.4p imply that it 7l = u 72 . 

(Hi) Finally, condition u 71 = u 12 is equivalent to Du^ = Du^ 2 , which holds 
as well if and only if 71 = 72 

First step is relatively easy since there is no regularity of 7 to preserve 
and thus one can extend by in D \ 0. Second step is dealt with in |11[ 
Proposition 6.1]. It is shown that H^ = 7i^ 2 implies /^(z, k) = f^ 2 (z,k) 
for all k € C and \z\ > 1. As a consequence (ii) follows. 

The step (Hi) is more complex because uniqueness results and a pri- 
ori estimates for pseudoanalytic equations in C like (|2.4|) only hold if the 
coefficients or the solutions decay fast enough at 00. Unfortunately the 
required decay properties for r seem to require roughly one derivative for 
7. However in [TT] it is shown that in the measurable setting at least we 
obtain subexponential decay. That is, we can write, 

Ul (z,k) = e **(*+«M(*.*)) (2.5) 

for some function e = e M (z, k) satisfying 

lim HfyfoAOHiwc) =0. 

This would not be enough if we would consider equation (|2.4|) for a single 
z. However, in [TT] it is used that u(z, k) solves an equation for each z. 
Further, one has asymptotic estimates for u both in the k (as above) and z 
variables. Then, a clever topological argument in both variables shows that, 
with these estimates, 77 determines the solution to (|2.4|) . 

Strategy for stability: In order to obtain stability, the natural idea is to 
try to quantify in an uniform way the arguments for uniqueness. This was 
done in [15] for C a conductivities. Let us recall the argument and specially 
the results which did not require regularity of 7 and would be instrumental 



for the current work. Let p = ||A 7l — A 72 ||. First one reduces to the unit 
disk by an argument which involves Whitney extension operator, the weak 
formulation (11.30 and a result of Brown about recovering continuous con- 
ductivities at the boundary (|19j). Next we investigate the relation between 
the corresponding scattering transforms. 

Theorem 2.3 (Stability of the scattering transforms). Let 71,72 be con- 
ductivities in B, with -^ < 7$ < K, and denote pi = jt^t- Then, for every 
k € C it holds that 

\r^(k)-r^(k)\<ce c ^p. (2.6) 

where the constant c depends only on K . 

The estimate is just pointwise but on the positive side it holds for L°° 
conductivities. In [151 Theorem 4.6] it is also given an explicit formula 
for the difference of scattering transforms which might be of independent 
interest. Next we state a result that is implicitely proved in [15, Theorem 
5.1]. There it is stated as a property of solutions to regular conductivities. 
However, in the proof the regularity is only used to obtain the decay in the 
k variable. Because of this, here we state it separately as condition (12.71) . 



Theorem 2.4 (A priori estimates in terms of scattering transform). Let 
K > 1 and 71, 72 be conductivities on D ; with -^ < 73 < K. Let 

1 j \ ik(z+e u . (z.k) ) 

u 7j (z,fe) = e *• *V '> , 

denote, as in (|2.5|) . the complex geometric optics solutions to ([Lip . Let us 
assume that there exist positive constants a, B such that for eack z,k £ C, 

Mz,k)\<^- (2-7) 

Then it follows that: 

A There exists new constants b = b(K), C = C(K,B), such that for every 
z € C there exists w € C satisfying: 



(a) \z-w\ <CB logi 

(b) u 7l (z,k) =u~f 2 (w,k). 



-ba 



where p = ||A 7l — A 



71 1 ^72 I 



B For each HC, there exists new constants b = b{K) and C = C(k, K) 
such that 

CBTi 

\\u~ fl (z,k) -u 12 (z,k)\\ L ^ {D4A{z)) < (2.8) 
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ba 


1 


K 


log- 




p 





Proof. The proof of A follows from |15|. Proposition 5.2] and [15} Proposition 
5.3]. Let us prove B. Given z G C, let w G C be given by part A. Then 

\u^(z, k) — u l2 (z, k)\ = \tiy 1 (z, k) — u 7l (iy, k)\. 

By the Holder continuity of -KT-quasiregular mappings, together with (a), we 
get 

\u yi (z,k) -u l2 {z,k)\ <C(k,K)\z-w\T< <C(k,K)CT< BT< 

and the desired estimate follows after renaming the constants. □ 

Unlike in the uniqueness arguments, going from u 71 — u 72 to D{u lr — u 12 ) is 
more delicated in the stability setting, since functions do not control their 
derivatives in general. This is solved in [T5], under Holder regularity, using 
the following fact. 

Theorem 2.5 (Schauder estimates). Let 7$, i = 1,2 be conductivities on D, 
such that -^ < 7j < K and |[7i||c a (B) < ^o- As always, denote fii = jr^r-, 
and let / w (z, k) be the corresponding complex geometric optics solutions to 
VTi\) . Then 

1. For each k € C there esists a constant C = C{k) > with 

\\U 1 (;k)-U 2 (;k)\\ cl+a{0) <C(k). (2.9) 

2. The jacobian determinant of /^(z, k) has a positive lower bound 

J(zj^(-,k))>c(K,k,r ). 

Now, to finish the proof of stability for Holder continuous conductivities, just 
note that an interpolation argument between L°° and C 1+a gives Lipschitz 
bounds for Df^. Thus, by \i = = and the second statement above, one 
obtains L°° stability for [i\ — [ii- The corresponding result for 71 — 72 comes 
due to fl2~2l). 



Strategy for stability under Sobolev regularity In the current work 
we will try to push the previous strategy to obtain I? stability. The previ- 
ous analysis shows that we can rely in many of the results from |114 [15] . In 
particular, we only have to prove that r M 1— > \i is continuous. 
For this, we start by reducing the problem in Section [3l We replace the 
assumption 7, G W a > p (Q) by 7, G W^' 2 (B), where < < min{i,a}. For 
this, it is used there that characteristic functions of Lipschitz domains be- 
long to W@' q (C) whenever f3q < 1. 
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Then we follow by investigating the regularity of solutions of Beltrami equa- 
tions with coefficients in fractional Sobolev spaces in order to obtain an 
estimate like (I2.9p . with the C l+a norm replaced by the sharp Sobolev norm 
attainable under our assumption on the Beltrami coefficient (see Theorem 
I4.6J) , It is also needed here to understand how composition with quasicon- 
formal mappings affects fractional Sobolev spaces. As far as we know, the 
estimates here are new and of their own interest. 

Afterwards we prove that our Sobolev assumption on fi suffices to get the 
uniform subexponential growth of the geometric optics solutions needed in 
condition (12. 7p in Theorem 12.41 (this is done in Section [5j see Theorem I5.7J) . 
In fact we obtain a very clean expression for the precise growth, achieving 
that the exponent depends linearly on a. Finally, in Section [6] we do the 
interpolation argument. Here we do not have enough regularity to control 
W 1,oc norms and here is where one sees why we need to be happy with the 
control on ||//i — A*2 1 1 z, 2 (K») ■ Also we do not have a pointwise lower bound for 
the corresponding Jacobians which causes also difficulties. 

3 Fractional Sobolev spaces and Reduction to /iG 

jy a ' 2 (B) 

3.1 On fractional Sobolev Spaces 

Following [1, p. 21], for any domain Q, we denote by W 1,P (Q,) the class of 
L p (£l) functions / with L p (f2) distributional derivatives of first order. This 
means that for any constant coefficients first order differential operator D 
there exists an L 2 (Q) function Df such that 



fD<p = - / Df<p 
n Jn 

whenever if G C°° is compactly supported inside of £1. Similarly one can 
define the Sobolev spaces W m ' p {Sl) of general integer order m > 1. 

It comes from the work of Calderon (see [21 p. 7] or [45]) that every 
Lipschitz domain Vt is an extension domain. That is, for any integer m > 
and any domain O' D O there exists a bounded linear extension operator 

E m : w p > 2 (n) -► w™' p (n') 

and therefore for every function / € W m,2 (Q) there is another function 
E m f 6 W m 'P(n>) such that E m f\ a = f. Of course, E m f G ^(C). 

Let us introduce for general domains $7 and any real number < a < 1 the 
complex interpolation space 

W a > p {ty = [17(0), W 1 '^)]^ 
12 



The closure of Cg°(f2) (C°° functions with compact support contained in 
0) in W a ' p (n) is denoted by W a ' p (fi). Functions in W a ' p (fi) can be 
extended by zero to the whole plane, and the extension belongs to W a ' p (C). 
Thus, we can identify any function in Wq (O) with its extension in W a ' p (C). 

When f2 is an extension domain, an interpolation argument shows 
(see [TJ p. 222]) that W a,p {Q) coincides with the space of restrictions 
to n of functions in W a ' p (C). That is, to each function u 6 W a > p (VL) 
one can associate a function u G W a,p (C) such that u^ = u and 

We have chosen just one way to introduce the fractional Sobolev spaces. 
In the rest of the subsection, we discuss the alternative characterizations 
and properties of these spaces needed in the rest of the paper. Two good 
sources for the basics of this theory are [U Chapter 7] , [lU Chapter 4] . 

Fourier side. For p = 2, it is easy to see that, 

W a >\C) = {/ G L 2 (C); (1 + |£| 2 )i /(0 G L 2 (C)} 

and that this agrees with the space of Bessel potentials 

W a >\C) = {f = G a *g;g£L 2 (C)} 

where G a is the Bessel Kernel [2j p. 10]. For p / 2 the situation is more 
complicated but it can be shown that 

W a *(C) = {/ G L p (C); ((1 + |£| 2 )i /(0) A G i p (C)| . 

Integral modulus of continuity We define the L p -difference of a func- 
tion / by 

w P (/)(y) = ||/(- + y)-/(-)llLp(o- (3-1) 

( see [35j Chapter V] Then the Besov spaces Ba q (W a ) are defined by 



B p f (R») = {/ G L P (M") : / Wp (/)(y) 



?| y |-(n+ag) < ^ 



There are many relations between Besov and fractional Sobolev spaces. We 
will need the following two facts, 

B% 2 = W a > 2 , W a ' p cB p ' 2 (p<2). (3.2) 

For a proof see [U Chapter 7] or [3SI Chapter V] . 
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Leibniz Rule [|31j] 

Lemma 3.1. Let a G (0, 1) and p G (1, oo). 
(a) Let f,g eCg°(C). Then, 

\\D a (fg) - fD a (g)-gD a (f)\\ p < C \\D^(f)\\ Pl \\D C K! , nu , 2 



whenever ax,ct2 G [0, a] are such that a± + ai = a and pi,P2 G (1, oo) 
satisfy 1 = -. 

Jd P\ PI p 

(b) Let /,^€Cg°(C). Then 

\\D a (fog)\\ p <C\\Df(g)\\ Pl \\D a g\\ P2 

whenever p 1 ,p 2 G (l,oo) satisfy i + i = |. 

fcj Let f,g£C^{C). Then 

\\D a (fg) -fD a (g) -gD a (f)\\ p < C \\D a (f)\\ p \\g\U 

whenever < a < 1 and 1 < p < oo. 

Remark 3.2. From property (a) and (c) it follows the generalized Leibniz 
rule 

\\D a (fg)\\ p < \\D a f\\ pi \\g\\ P2 + \\D a g\\ P3 \\f\\ P4 (3.3) 

whenever 1 < Pi,P2,Pr,P4 < oo and - = 1 = 1 . Moreover if 

— rm-iiroiri — p pi P2 PA Pi 

suppt(/, g) G D we have that 

\\D a (fg)\\ p < \\D a f\\ Pl \\g\\ P2 + \\D a g\\ LP3(p) \\f\\ L v 4(D) (3.4) 

Pointwise Inequalities 

Lemma 3.3 (Pointwise inequalities, p]). If f G W a ' p (C), a > 0, 1 < 
p < oo, then for each < A < a there exists a function g = g\ G L Px (C), 

P* = 2-(l-X)p SUch that 

\f(z) - /HI <\z- w\ x (g(z) + g(w)) (3.5) 

for almost every z,w G C. Furthermore, we have that 

IIsIIlpa(c) - CaII/IIw«.j'(c)> 

and i/ie constant C\ remains bounded as A — > a. 
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3.2 Reduction to p = 2 

This reduction relies on the fact that \i € L oc (C)nW a ' p (C) and the following 
interpolation Lemma. 

Lemma 3.4. Let f e W a °' po nW ai ' Pl , where 1 < p ,p± < oo, < a ,«i < 
1, and 9 G (0, 1). Then, 



< 



where 

1 O 1 fl 

a = 9 oto + (1 — 9) a\ and — = 1 . 

P Po Pi 

Furthermore, if either p$ = oo or p\ = oo, then the above inequality holds 

true by replacing W ai ' Pi by the Riesz potentials space I ai * BMO. 

Proof. It is well known that the complex interpolation method gives 

whenever 1 < p < oo (for the proof of this, see for instance [51]). For 
p = oo, the same result holds true if we replace W a,oc by the space of Riesz 
potentials I a * BMO of BMO functions (for this, see [33]). □ 

Let /ibea compactly supported Beltrami coefficient. Then, it belongs both 
to L l {C) and L°°(C). If we also assume that /i € W a ' p (C) for some a,p, 
then we can use the above interpolation to see that fi 6 W^' 9 (C), for any 
1 < q < oo and some < (3 < a. We are particularly interested in q = 2. 

Lemma 3.5. Suppose that \i € W a,p {Q) n L°°(yt) for some p > 1 and 
< a < 1. Then, 

• For any < 9 < 1, 

1 1 1 1 , ,, ,,1 _/i . . . . a 

IMIw^.f (n) ~ IMIi,<»(n) IMIw^n)' 

• For any < 9 < 1, 

.... .. ..1 a . . . . a 

• One always has 

n n ^ m(Ts \ ii IIP*/ 2 

where (3 = ^- and p* = min{p, ^j}- 

Proof. The first inequality comes easily interpolating between BMO(£l) and 
W Q ' p (fi) (see |43] for more details). For the second, simply notice that 
compactly supported Beltrami coefficients belong to all L P (Q) spaces, p > 1, 
so one can do the same between L 1+£ (f2) (e as small as desired) and W a ' p {Vl). 
The last statement is obtained by letting 9 = \ above. □ 
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3.3 Reduction to Q = ED and \i e Wg' p {B) 

The proof of the following lemma relies in the fact that characteristic func- 
tions of Lipschitz belong to W a ' 2 for each a < ^. 

Theorem 3.6. Let tt be a Lipschitz domain, strictly included in B. Let 
H G W a ' 2 (n). Define 

_ _ (V n 
^ ~\o c\n" 

Then, ft G Wq (C) /or /3 < min{a, -} and 

llAllw/3>2(c) < C|IHIw a ' 2 (C)- 

Analogous results can be stated for the extensions by 1 o/7j. 

Proof. Since O is an extension domain, there is an extension /xo of [i belong- 
ing to W a ' 2 (C). Of course, such extension //o need not be supported in Q 
any more. Now \i can be introduced as the pointwise multiplication 

By virtue Lemma 13.11 it is enough to study the smoothness of the character- 
istic function xu- A wa Y to see this is to recall that fraccional Sobolev spaces 
are invariant under composition with bilipschitz maps |54| . Now, the char- 



acteristic function of the half plane belongs to W^f(C) whenever ap < 1. 
Therefore, by a partition of unity argument, we get that xn € W a,p (C) 
when ap < 1. The proof is conclude. □ 

Now we need to compare the original Dirichlet-to-Neumann maps with the 
Dirichlet-to-Neumann maps of the extensions. 

Lemma 3.7. Let Q be a domain strictly included in D. Let 71,72 G L°°(Q) 
be conductivities in fl. Further, assume that 

^ < li{z) < K 

for almost every z G £1. Let ji denote the corresponding extensions by 1 to 
all of C. Then, 



I 71 72 II u^ fan\^,tj-^ /an\ — II 71 72 



1 



Proof. We follow the ideas of |X5|, Theorem 6.2], although the stability result 
from [TH] is not needed in our situation. Let ipo G H^(dB>). Let Uj G -ff 1 (D) 
be the solution to 

J V • (TjVu,) = inB 

1 Uj = ipo in SB. 
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Let also «2 be defined by 



I V ■ (72VM2) = in O 
I ui = U\ in dQ. 

Define now V2 = 1*2 Xfi + u\ Xo\ti- As in [15], we first control U2 — v<i in terms 
of p. To do this, 



V(v 2 ~U2)\ 2 <C I 72 V(u 2 - U2) • V(t>2 - «2) 



c / 72 Vv 2 • V(v 2 - U2) 



because V2 — U2 £ i?d(B) and the 72-harmonicity of U2 in B. By adding and 
substracting J^tlVwi • V(ii2 — U2), and using that 71 = 72 = 1 off O, the 
right hand side above is bounded by a constant times 

/ 7iVtii • V(i)2 — M2) + / (71 Vui - 72 Vu 2 ) • V(u2 - "2) 

Here the first term vanishes because u\ is 71-harmonic on B and £2 — U2 G 
i/g(B). For the second, we observe that u\ is 71-harmonic in f2, U2 is 72- 
harmonic in Q, and 1*2 — «i € i^g(O). Thus, 



/ (71 Vui - 72 Vu 2 ) • V(v 2 - u 2 ) 
Jn 



= |((A 7l - A j2 )(u lldn ),(v 2 - U2)\an)\ 

< p lluiH 1 \\vo — Uoll 1 

<p||Vni|| L2(n) ||V(w2-n 2 )|| L 2 {Q) 



Summarizing, we get 

|V(«2 --U 2 )| 2 ) < Cp || Vttl ||i2(0) < cp||Vui|| L 2 ( 
< cp Hvoll ! 



(3.6) 



We will use this to compare the Dirichlet-to-Neumann maps at <9B. If ipQ S 
H^{dJ}) is any testing function, and ip is any i7 1 (B) extension, 



(3.7) 



((A 7l - A 72 )(<^ ), V>o) = / (71 V«i - 72 Vn 2 ) • V^ 



We will divide the bound of this quantity in two steps. For the first, 



(71 Vui - (72 Xn + 7i Xn\n) ^v 2 ) ■ Vip 



|((A 71 -A 72 )(w 1 | an ),V|9n>| 
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which is bounded by 



p||«i|9nll K i (an) ll^n|l H i (9n) <p||Vfii|U2 (n) ||V^|| L 2 {n) 



< pI|Vmi|| L 2 (d) ||V^||z,2 ( 



< P Vo l 



,J0 " H h 



We are left with, 



/ ((72 Xn + 71 Xb\q)V£ 2 - 72 Vu 2 ) • VV> 

JO 



which is smaller than, 



72 V(w 2 - M2) • V-0 + / V(£ 2 - -u 2 ) • V^ 

<> Jo\n 



which in turn is controlled, using (|3.6p . by a multiple of 

|V(«2 " «2)| IW>| < ||V(V2 " U 2 )|k2(D) ||V^|| L 2 ( 



/ 



< CO I to I 1 

— r ll ^ ll H2(aiD)) 



^0 1 
M Jf^(9]D)) 



This gives for (|3.7p that the difference of Dirichlet-to-Neumann maps satis- 
fies 

l(( A 7l - A 7 2 )(^0),V'0)| <Cp\\<p\ 



<H1(8B) 



U|l _ff^(eiD)) 



as desired. 



D 



Remark 3.8. The trivial extension of the conductivities by 1 simplifies the 
arguments but has the prizes or loosing regularity if a > 1/2. An argument 
similar to that in |15j would need an 1? version of the boundary recovery 
result of Brown (see also [5]) of the type 

II71 - 72||L2(an) < Cp 



4 Beltrami equations and fractional Sobolev 
spaces 

This section is devoted to investigate how quasiconformal mappings inter- 
play with fractional Sobolev spaces. We face three different goals. First, 
given a Beltrami coefficient \x € Wq' (C), we find j3 G (0, a) such that for 
any X-quasiconformal mapping <j> the composition fj, o (f>, which is another 
Beltrami coefficient with the same ellipticity bound, belongs to W^' 2 (C). 
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Secondly, we obtain the optimal (at least when a « 1), Sobolev regularity 
for the homeomorphic solutions to the equation 

under the assumptions of ellipticity and Sobolev regularity for the coeffi- 
cients. Finally, we obtain bounds for the complex geometric optics solutions. 
Many properties of planar quasiconformal mappings rely on two precise in- 
tegral operators, the Cauchy transform, 

and the Beurling transform, 

Tip(z) = — lim / . ^ W \, dA{w). (4.2) 

Their basic properties are well known and can be found in any reference 
concerning planar quasiconformal mappings, [3l \9\ \TT\. 

4.1 Composition with quasiconformal mappings 

Let n be a compactly supported Beltrami coefficient, satisfying 

K-\ 

Further, assume that 

H 6 W a ' 2 (C) and ||A*||wa,a ( Q < T 

for some a > and some To > 0. Let eft : C — > C be a planar K- 
quasiconformal mapping. In this section, we look for those /3 > such that 

flO(f>£WP' 2 {C). 

We need to recall a local version of a lemma due to Fefferman and Stein, 
see |40| and |26|. Proposition 2.24]. The proof follows from Vitali covering 
Lemma, exactly as in |40| . By Mf we denote the Hardy-Littlewood maximal 
function, 

M/(x) = sup-^ / /. 

where the supremum runs over all disks D with x G D, while Mcif denote 
its local version, that is, 

M n /0*0 = sup-t- // 

l-^l JD 

where the supremum is taken over all discs D with x € D C fi. 

19 



Lemma 4.1. Let w > a locally integrable function. Then 

[ \M n f\ p iodx < [ \f\ p Mu. 

We can now prove the main result of this section. 

Proposition 4.2. Let K > 1. Let fx £ W a ' 2 (C) for some a € (0,1), 
and assume that \fi\ < -^tjXid- Let 4> '■ C — > C be any K-quasiconformal 
mapping, conformal out of a compact set, and normalized so that \(j)(z)—z\ — ► 

as \z\ — ► oo. T/ien 

whenever (3 < j^ . Moreover, 

j_ 

||/X O ^>|| W /3,2( C) < C |M|^ a , 2(C) , 

/or some constant C > depending only on a, (3 and K. 



Proof. It is clear that /i o <^> belongs to L 2 (C), so since W a ' 2 agrees with the 

ci 



Besov space B a ' , it suffices to show the convergence of the integral 



IcJc H 2+2/3 



dA{z) dA(w) 



for every j3 < ^. First of all, note that by Koebe's j Theorem we have 
the inclusions </>(D) C 4B and c/> _1 (4B) C 16D. Thus, for large w there is 
nothing to say since 



J r\ji^{z + w ))-n{<t>{z))\ 2 

J\w\>l JC 



.., ,, , dA(z)dA(w) 



<2||M||| 2(C) / 1 ^^H = ^f !e 
j \w\>\ \ uj \ y 

for some universal constant C > 0. Then we are left to bound the integral 

/^, + .)M»))f 

M<iJc M 2 + 2 ^ 

Notice that this integral is in fact the same as 

/ Q2T2^ dA(z)dA(w) (4.3) 

|tu|<l JF \ w \ H 

where F = {z € C : d(z, $> _1 (D)) < 1} C 17B by Koebe's Theorem again. 
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Since y, G W a ' 2 (C), then by interpolation we get that fj, G VF ae, e(C) 
for each G (0, 1) and 

W^w ae -i{C) -^Nlw IMIw<*. 2 (C)- 

The goal is to choose to obtain larger possible /?. For this, we use condition 
(I3.5p , Indeed, for each A G (0, aff) there exists a function g = g\ G L PA (C), 
Pa = g„ (a 2 g_ A) , such that 

MO ~ m\ < \C ~ Z\ X (9(() + 9(0) 
at almost every (,(GC. Furthermore, 

<r r> ll,/l|i-s ll/,ll e 

-i °A HMlloo HMllvK a . 2 (c)- 

with Ca bounded as A — > a#. Hence, if \w\ < 1 then 

W(* + «0)-Mfl*))l < p(^)-<^)l\ A W(z + w)) + 5(0(z))) . 

Now use quasiconformality and the reverse Holder inequality for the jacobian 
(see [10] for a more precise formulation) to get that 

(z + w) - 4>(z)\\ / diam (f>(D(z, \w\))^ ' 



w\ / \ diamZ?(z, |u)j 



— /" 



< c k rT^7 i im / J(C,#dA(C) 



l-DU, I'WDI ./£>(*,|io|) 



< C K (M n J A ^)) 5 

where = {z6C:d(x,f 1 (D))<2} and M n J x (z) denotes the local 
Hardy-Littlewood maximal function Mq at the point z of J(-,<p) x . Note 
also that f2 C 18B by Koebe's Theorem. By symmetry, we could also write 
MqJ\(z + w) instead of MqJ\(z), so we end up getting 

\lM(<f>(z + w))-fi(tf>(z))\ 2 



III.' 



2A 



< C {M n J x (z + w) g(<f>(z + w)f + MnJx(z) g^(z)f) . 



Therefore the integral at (|4.3|) is bounded by a constant times 
r r M n J x (z + w)g (<t>(z + w)) 2 + M n J x (z) g((t>(z)) 2 

J\w\<\ JF 



w |2+2/3-2A dA ( Z ) dA ( W ) 



If we restrict ourselves to values of A within the interval (f3,a0), then the 
integral above is bounded by 

C 



, M ft J A (z)^(^))^dA(z). (4.4) 
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To get bounds for this, we start by choosing parameters. Fix a, (3 and K 
with (3 < -^ . Then we can find s > 1 such that (3 s < j^ . Now let us consider 
real numbers 6 £ (0, 1) and A £ (/3, a9) satisfying 

A + (l- a )0<-L (4.5) 



Such conditions are compatible precisely when (3s < #. Condition (|4.5p 
also guarantees that p\ > 2Ks, so we can find r satisfying 

l + Xs(K-l)<r<^(l + Xs(K-l)), (4.6) 

IK s 

and obviously r > 1. 

Once the parameters have been chosen, we proceed as follows. First, by 

Holder's inequality 



M n Jx(z) g o <P(zY dA(z) = / M Q Jx(z) g o <p(zy xf(z) dA(z) 

< (J^M n J x (z) s go^zf s XF(z)dA{z)y |F|H 
Now we use Lemma [4. II to bound the last integral above by a constant times 

J(z,<A) As Af( 5 o^ XF )(z)^(z). 



in 
For any r > 1, we can bound the above integral by 

i . - 



J(z, 0) AS M(g o «^ XF )(zf cM(z) / J(z, <p)^dA(z) 
c / \Jn 

The first inequality at (14. 6p guarantees that the weight J(-,0) belongs to 
the Muckenhoupt class j4 r (see [TU] for details). Therefore we can use the 
weighted L r inequality for the maximal function and a change of coordinates 
to see that 

\As -\r ( i2s „, \ t „\r j \ i „\ ^- r< I it _ j.\As „ „ j./~\2sr 



J(z, 0) As M( 5 o «^ s XF )(z) r <L4(z) < C r y J(z, 0) AS 9 o <t>( z y sr dA(z) 

= C r f J(w,(f)- l ) l ~ Xs g(w) 2sr dA(w), 

U{F) 

where C r is a positive constant depending on r. Summarizing, we get for 
the integral at (14. 4|) the bound 

C|F|H; (f J(z,cf>) Xs dA(z)) a ' if J(w,r l ) l - Xs g(w) 2sr dA(w)Y 
\Jn J \J<t>(F) J 
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For the second integral above, we only use Holder's inequality again, and 
obtain the bound 



2r.- 



g(w) Px dA(w)\ / J(w,<t>~) p *~ 2rs dA(w)\ 

0(F) J \JcKF) J 

which is finite provided that both p\ > 2rs and Vx _ 2 rs ^ Tf^T hold. But 
both facts are guaranteed by our choice of parameters, in particular to the 
second inequality at (|4.6H . This means that the integral at (|4.4p has the 
upper bound 

^ ^J^{zA) Xs dA{z)^ lbll!pA W F)) ( / ^ A^ -1 ) p ^ s rfAMj 

Since both ^ and _1 are normalized i^-quasiconformal mappings, the two 
integrals above are bounded by constants depending only on K. One obtains 
for the integral at (|4.3p the bound 



LAI 



C C 



l- 



— — tII5Hlpa(0(F)) < — — — T IIMIIi«(c) \W\w«fl^j 

where the constant C depends on r, s, A, 9, a and K. 

To find larger possible /?, we have to find the supremum of those A for which 

the pair (8, A) belongs to the set 

A = 1(9, A); < 9 < 1, /5 < A < 9a, \+(l-a)9<-^- 
{ Ks 

according to (|4.5|) , This supremum is easily seen to be -^-. Further, all the 
above argument works for every s 6 (1, -M?), so that the bound for (|4.3p 
reads now as 

H 1 



K r 
Summarizing, 



// 






C JC 

/II//II 2 1 

~ M /3 + a -/3 l^ll^,2 (€ ) 



2 



for all /? G (0, #). Equivalently, we have an inequality for the nonhomoge- 



K 

neous norms 



\\H o 0|| W /3,2( C ) < C(a,/3,K) ||HI^,2 (C) , 
as stated. D 
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Remark 4.3. The condition (3 < % is by no means sharp. This is clear 
when q is close to 1 but also it can be seen from the fact that we are using 
the Holder regularity of (ft. It seems that one can get a factor which runs 
between 1/K and 1 in terms of a. As promised in the introduction this will 
be a matter of a forthcoming work. 

4.2 Regularity of homeomorphic solutions 

We start by recalling the basic result on the existence of homeomorphic 
solutions to Beltrami type equations. In absence of extra regularity the 
integrability of the solutions comes from the work of Astala [7J. We recall 
the proof in terms of Neumann series since it will be used both in this section 
and in the sequel. 

Lemma 4.4. Let /i, v be bounded functions, compactly supported in D, such 
that \\fi(z)\ + |^(z)|| < -j^ri at almost every z € C. The equation 

df = fidf + udj (4.7) 



admits only one homeomorphic solution (ft : C — > C, such that \<ft{z) 

2K - 
K-\. 



0(\j\z\) as \z\ — ► oo. Further, if p £ (2, K _ 1 ) then the quantity 



\\9(J>-1\\lp{C) + [|00||lp(C) 
is bounded by a constant C = C(K,p) that depends only on K and p. 
Proof. Put (ft{z) = z + Ch(z), where h is defined by 

(I - [iT - vT)h = fi + u. 

and C and T denote, respectively, Cauchy and Beurling transforms. Since 
T is an isometry in L 2 (C), one can construct such a function h as Neumann 
series 



/i = X>T + I /T) n ( M + 



n=0 
It 



which obviously defines an L (C) function. By Riesz-Thorin interpolation 
theorem, 



lim \\T\\„ = 1, 

p^2 



K+l 



it then follows that h £ L P (C) for every p > 2 such that ||T|| P < 
Hence, the Cauchy transform Ch is Holder continuous (with exponent 1 — ). 
Further, since h is compactly supported, we get \(ft(z) — z\ = \Ch(z)\ < l^j, 
and in fact (ft — z belongs to W 1,P (C) for such values of p. A usual topological 
argument proves that is a homeomorphism. For the uniqueness, note that 
if we are given two solutions (fti, <p2 as in the statement then d((ft\ o(ft^ ) = 
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so that <pi o 2 (z) — z is holomorphic on C and vanishes at infinity. 

In order to recover the precise range (^xj, -j^zj) obtained by Astala [7] we 

need a remarkable result from [TO] which says that I — fiT — vT : L P (C) — » 

LP(C) defines a bounded invertible operator for these values of p. Further, 

both its LP norm and that of its inverse depend only on K and p. This 

implies that for every p 6 (^xjj ^zj) there is a constant C = C(K,p) such 

that 

||/i||jy(C) < CW, p . 

The claim follows since d(ft — 1 = Th and dcf) = h. □ 

Once we know about the existence of homeomorphic solutions, it is time to 
check their regularity when the coefficients belong to some fractional Sobolev 
space. 

Theorem 4.5. Let a € (0, 1), and suppose that fj,, v € W a ' 2 (C) are Beltrami 
coefficients, compactly supported in D, such that 

IMOI + KOII < §^- 

at almost every z € B. Let <j) : C — ► C be the only homeomorphism satisfying 

d(j> = fidcj) + vdcj) 

and 4>(z) — z = 0(1/ z) as \z\ — » oo. Then, 4>(z) — z belongs to W l+Sa,2 (C) 
for every 9 E (0, -g), and 

\\D 1+ea (cf) - z)\\ L 2(C) < c k (||mI!h^,2( C ) + M\ e wa ,2(c)j 
for some constant Ck depending only on K. 

Proof. We consider a C°° function ijj, compactly supported inside of D, such 
that < tp < 1 and ftp = 1. For n = 1,2,... let tp n (z) = n 2 ip(nz). Put 



and 



Hn(z)= \ H(w)lp n (z -w)dA(w), 

Jc 



v n (z)= \ v(w)i) n (z-w)dA(w). 
'C 



It is clear that both /i n , v n are compactly supported in ^^-U, |Mn( z )l + 

K(^)| < ^1, WfJ-n ~ H\\w a ' 2 (C) -> and IK ~ u \\w a < 2 (C) "^ as TH OO. 

Indeed there is convergence in LP for all p € (1, oo). Thus, by interpolation 
we then get that for any < < 1 

lim \\u n — /ill „„ 2 + IIia, — i/|| „ fl 2 =0 
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and in particular, the sequences D \i n and D v n are bounded in L«(C). 
Let (j> n be the only if-quasiconformal mapping <p n : C — ► C satisfying 



04> n = Vn d<j> n + v n d(j) n (4.8) 

and normalized by 4> n ( z ) — z = O n (l/z) as \z\ — ► oo. By the construction in 
Lemma |4.4( <j) n (z) = z + Ch n (z) where /i n is the only L 2 (C) solution to 



/l n = fj, n Th n + U n Th n + (jJLn + v n ), 

and Ch n denotes the Cauchy transform. As in Lemma 14.41 h n belongs to 
L P (C) for all p € (^xj, ^xi) an d ||/in||z,J>(C) < C = C(K,p); in particular, 
n — 2 is a bounded sequence in PF > P (C). 



We now write equation (|4.8p as 



<9(<^ n - z) = H n d{(/)n - z) + ^„<9(> n - z) + jU„ + u n 

and take fractional derivatives. If (3 = a6, we can use Lemma l3.ll (a) to find 
two functions Ep, Fp such that 

D^n - Z) 

= Z)V d{cf> n -z)+fi n D p d{cf> n -z)+Ep 

+ d V a(^„ - z) + i/ n £>^a(0 n - z) + f^. 

Further, E^ satisfies 

WEpWviC) < C \\D^\\ LV1{C ) \\d((f>n ~ z)\\ LP2{c) , (4.9) 

where p2 is any real number with 2 < p2 < -^i an d ~ + = 2' an< ^ Co 
depends on pi,p2- Analogously, 

ll^ll^(C) < Co P^lk^c) ll^n " 2)IIlw(C)- (4-10) 

Now we notice that we have D^dip = dD^cp and similarly for 9. Further, if 
<£> is real then D^ip is also real. Thus 

dD^(4> n - z) 
= n n dD^{(j> n -z)+ D^ n d{(j> n -z)+E p 



+ v n d{B^ n - z)) + D?v n d^n -z) + Fp, 
or equivalently 

(/ - fi n T - u n T)(d) (DP&n - z)) 

= L> V 9{(j) n -z) + D?v n d((j) n - z) +Ep + Fp. 

26 



The term on the right hand side is actually an L 2 (C) function. To see this, 
it suffices to choose in both fj4.9f) and (|4.10p the value pi = 1 for some 
6 € (0, ^). Now, the operator I — /j n T — u n T is continuously invertible in 
L 2 (C), and a Neumann series argument shows that the norm of its inverse 
is bounded by ^(K + 1). Thus, 

pi^(0 n -z)|| L2(c) 

i^r + 1 /K-i\ 1_e n \ 

< ^0 2 I K + 1 / (JI*Mw a - 2 (C) + ll^n lw«< a (C) ) \\°K < Pn—Z)\\LP2(C) 

< Co ^ (^llMn|lvy«.2(c) + ll^nll^a.afoj W^i^n ~ z)||lp2(Q 

where Co is the constant in (|4.9p . As n — > oo, the right hand side is bounded 
by 

Co — 2 — (HmHw«. 2 (C) + IHw«- 2 (C)J C (K,p 2 ) 

because \\d(4> — z)\\ LV 2(C) < C(K,p2). Hence, dD"((j) n — z) is bounded in 
L 2 (C), and thus also dD@((j) n — z), because T is an isometry of L 2 (C) and 
T(dD f3 ((f> n — z)) = dD f3 ((f> n — z). Therefore, by passing to a subsequence 
we see that D^{<j) n — z) converges in W 1,2 (C), and as a consequence <fi — z 
belongs to W 1+ ^ ,2 (C). Further, we have the bounds 

II nl+Sa/i _A|| s /~i Ml, ,110 i II, ,110 I 

\\ D 10- z )\\l\C) < ^ ywWw^iC) + IMIiy«> 2 (C) ) 

for some constant C depending only on K. □ 

4.3 Regularity of complex geometric optics solutions 

We are now ready to give precise bounds on the Sobolev regularity of the 
complex geometric optics solutions to the equation df = fidf introduced in 
Theorem E2 

Theorem 4.6. Let p, € W a,2 (C) be a Beltrami coefficient, compactly sup- 
ported in D, with \\fi\\ oo < -fex and \\^\\w a > 2 {C) ^ ^V Let f = f^(z,k) the 
complex geometric optics solutions to the equation 

For any < 8 < -j* we have that 

feW? + c ea > 2 (C). 
Further, we have the estimate 

\\D 1+ae (/„)(• , k)\\ LHB) < C{K) e°\ k \ (1 + \k\) (To + \k\ a ) 9 
were C,C(K) > 0, and C(K) depends only on K. 
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Proof. The existence and uniqueness of the complex geometric optics solu- 
tions comes from [11| Theorem 4.2] (see Theorem 12.21 in the present paper). 
Secondly, it is shown in |114 Lemma 7.1] that / may be represented as 

f = e ik<i> 
where eft : C — ► C is the only W,' (C) homeomorphism solving 

d0 = -fi^e^ k (0)d0 (4.11) 

and normalized by the condition \cp(z) — z\ — > as \z\ — > oo. Here 
e_fc(0(z)) = e -* fc 0(*O-* fc 0(*) i s a unimodular function so that \e-k(<p(z))\ = 1. 
In particular, is conformal outside of B, because supp(/u) C B. Thus by 
Koebe \ theorem 

0(B)c4B=H|0(-,fc)IU«»(D)<4 (4.12) 

Our first task is to determine the Sobolev regularity of the coefficient 
//e_fc(</>) in equation (|4.11|) . We will argue by interpolation. Firstly, by 

||e- fc (0)|| L2(B) < |B|i 
For the I? norm of the derivative, we invoke Lemma 14.41 to obtain that 

\\D^-z)\\ L2{c) <C(K), 

Thus, 



\\D(e- k (<f>))\\f<p)<C(K)\k\ 
and by interpolation we arrive to, 



i ,.. a 



\\e- k (<l>)\\ w «* m <C(K)\B\2\k\ 

Now we will use the remark [3.21 to see that that /ue_fc(0) belongs also to 
W a ' 2 (C). Since e-k(4>) is unimodular, the L 2 bound is obvious. By virtue 
of (|3.4|) we have that 

||Zn/ze_ fc (0))|| L2(c) < C(\\D Q v\\ L2{c) \\e„ k (0)\\ + Kp Q ( e _ fe (0))|| L 2 (D) ) 

<C|B|^(|fc| a + r ) 

(4.13) 
The bound (|4.13p allows us to apply Theorem 14.51 to equation (|4.11[) . We 
obtain that <t>o(z) = <p{z) — z satisfies the estimate 

\\D ae (dfo)\\ L 2 {c) + \\D a9 (dfo)\\ LHc) < C K (To + \k\ a ) 9 (4.14) 

for € (0, j^). We push this bound to /. Since f(z) = e %k ^ z \ we have 

df( z ) = e ik<t) ^ ikd<p(z) 
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and again from Lemma 13, II (c), for any disk D, 

\\D a0 (df)\\ L2{D) <\\D ae {e ik ^ikd^\\ L , {D) + \\ikD* e (d<t>)e ik *\\ LHD) 

+ C\\e il *\\ LOO(D) \\D a0 (d<f>)\\ I ? (D) . * ■"' )) 

For the second and third terms on the right hand side above, we notice that 

(ETT21) yields that 

sup|e* fc< ^)| <e 4|fc| , (4.16) 

zee 

which combined with (|4.14p provides us with the estimate 



likD^id&e^h^ < \k\e°\ k \ \\D ae (dd>)\\ L2{ 

<\k\e c ^C K (T + \k\ a 



(4.17) 



Concerning the first term in (|4.15p . we recall that LP bounds for d(j> are 
possible only for p G C^rp ~^j)- This forces us to look for L q bounds for 
D a9 {e tk< f > ), for some q > 2K. These bounds are easily obtained by interpo- 
lation. More precisely, we know the L°° bound given at (|4.12p . Further, we 
have also a W 1 ' 2 bound, 

\d(e ih K*))ft dA(z) < \k\ 2 e c ^K [ J(z,<f>)dA(z) < \k\ 2 e c|fc| K |D|. 

Jo 

Thus, by interpolation we obtain 

II r\aO / ik4>\\\ <s r 1 \\Jk(j)\\l— a.6 \\ a( ik<j>\ \\a6 

i|jD (C }II L^(B)- C||e Hl«(D) H S ( e ^)lb(B) 

< \k\ ae e°\ k \ {K\B\)f < C(K) \k\ ae e c ^. 

Now recall that < 8 < -g is fixed, and let p = \. If we now consider any 
real number s such that l 2 ae < s < -J^r, we obtain 

n^(e^) ikm»<n < 1*1 n^(^)ii L - (D) m\ LT ^ m 

= C{K)\k\ 1+ae e c \ k \ 

because the normalization on </> forces uniform bounds for ||9^>||ls(]D)) depend- 
ing only on K. Summarizing, (|4.15|) gives us the bound 

\\D ae {df)\\ L 2 m < C{K) e c ^ (1 + |A;|) (T + \k\ a ) e . 

Similar calculations give the corresponding bound for D a9 (df). □ 

We will also need the following bounds in Section [6j 
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Lemma 4.7. Let fi be a Beltrami coefficient, compactly supported in D. 
Assume that ||/x||oo < kTi an ^ llMllw a ' 2 (C) — ^o- Let f = fn(z,k) denote 
the complex geometric optics solutions to 

df = »df. 

Let p < 2/(K — 1). Then, for any disk D 

Jd 



p 
<C 



df(z) 
where the constant C depends on diam(D) ; k, K and Tq. 

Proof. The function / can be represented as 

t = e ik<t> 

so that 

1 1 1 



df e ik <t> ikd<j>' 

As we have seen in (|4.16p in the proof of the above Theorem gives also lower 
local uniform bounds for e $ . Thus, only LP bounds for dcj) are needed. But 
these bounds come from the fact that is a normalized K-quasiconformal 
mapping, so that 

dc/)eL p 

□ 

5 Uniform subexponential decay 

We investigate the decay property of complex geometric optic solutions to 
the equation 

dfx = X^dfx, 

where A G <91D> is a fixed complex parameter, and /i G Wq' (C) is a Bel- 
trami coefficient compactly supported in B. It turns out that f\ admits the 
representation 

f x (z,k) = e ikMz ' k) 

where 4>\ satisfies the following properties (see |lll Lemma 7.1] or the proof 
of Theorem 14.61 above) : 

1. 4>\(- ,k) : C — > C is a quasiconformal mapping. 

2. (p\(z, k) = z + Ok(l/z) as \z\ — ► oo 

3. <j)\ satisfies the nonlinear equation 



d(f) X (z,k) = -X fj,(z) -e^ k ((f> x (z,k)) d(/>x(z,k) (5.1) 
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|&(*,*)-*|< nnE ( 5 - 2 ) 



As was explained in Section 2, our goal is to obtain a uniform decay of the 

type 

C_ 

\W 

The precise statement can be found at Theorem 15,71 For the proof, we 
will mainly follow the lines of both [111 [T5] . This consists on investigating 
first the behaviour of linear Beltrami equations with the rapidly oscillating 
coefficients fj,(z) e-k(z), and then treat the nonlinearity as a perturbation. 

5.1 Estimates for the linear equation 

As usually, \x denotes a Beltrami coefficient, compactly supported in B, with 

the ellipticity bound 

K-l 
llMllL-(n) < J^r[- K 

and the smoothness assumption 

IImIIw»' 3 (C) = IHIl 2 (c) + H^VIIi^C) ^ r o 

for some < a < 1 and Tq > 0. For each complex numbers k G C and 
A 6 B, let ip = ip\(z, k) be the only homeomorphic solution to the problem, 



(dip(z,k) = %\e_ k (z)n(z)dip(z,k) 
\ip(z,k)-z = 0(l/z), z^oo 

Then ip can be represented by means of a Cauchy transform 



(5.3) 



tp(z,k)-z = / dip(w,k)$(z,w)dA(w), (5.4) 

•Ac 

where $>(z, w) = z-w ^ or a smo °th cutoff function ipo = 1 on B (in partic- 
ular on the support of dip). We need subtle properties for both terms. 

Lemma 5.1. Let uq be given, and let s > 2 be such that 

k\\T\\ 3 < 1. 

There exists a decomposition dip\(z,k) = g\(z,k) + h\(z,k) satisfying the 
following properties: 

1. \\hx(;k)\\ s <C(K,s)(K\\T\\ s ) no . 

2. ||<7A(',fc)||. <C(k). 
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3. For R>0 and \k\ > 2R, 



" r 



\gx(£,k)\ q dA(£) < C(a,K,p) M(p) n ° -— 
£\<R / II 

where 1 < p < 2, q = ^ =l , g\(£,k) = {g\(-,k)) A (£) and the value of 
M(p) is given in 115. 8\) . 

The proof will rely on the Neumann series expression of dip. For this, we 
consider the unimodular factors 

e n (z) =e i<kz+ikz)_ 

An idea which goes back to [TT] is dealing with the unimodular factors e n 
conjugating them with the Beurling transform. Namely, we express 

(k \" +1 
d ^ = Yl ( k X ) e -(«+i) /*> ( 5 - 5 ) 

where 

[ f0 = fl (5.6) 

[fn = M r n(/n-l), n = 1, 2, . . . 

Here by T n we denote a singular integral operator defined by the rule 

T n (ip) = e n T{e- n (p) 



where T is the usual Beurling transform (|4.2p . It is not hard to see that T n 
is represented, at the frequency side, by a unimodular multiplier of the form 

§ — n 



Thus, 



l r n||L 2 (C) = l!^n||i 2 (C)^L 2 (C) 



and T n is an isometry of L 2 (C). In fact, for any 1 < p < oo, 

11^(^)11^(0) = \\T( e -n¥)\\LP(C) < \\ T \\lp(C) IMIi>(C) 

because |e n (z)| = 1, so that ||T n || iP ( C ) = ||T|| iP ( C ). As T n is given by a 
Fourier multiplier, it commutes with any constant coefficients differential 
operator D and thus, 

ll^n^llw/i'^C) = ll^n¥>||LP(C) + \\T n (D<p) IIlp(C) < ll^llp IMIw^C) 

and therefore ||T n ||^i, P (c) ^ II^IIlp(C)- Furthermore, the complex interpola- 
tion method gives that for any < j3 < 1, 

||2n.||yi//3>P(C) < Co II^IIlp(C) (5-7) 
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where Co > is a universal constant. 
Let 1 < p < 2 be fixed. We declare 



B = B{p) = ||T n ||wa,p(c) < \\T\\lp(q 
D = D(p) = \\T n \\j 

M = M(p) = B + D 



D = D{p) = \\T n \\ 2p <\\T\\ 2p ( 58 ) 



It is well known that the series Yl n fn defines a compactly supported L 2 (C) 
function (actually L P (C) for any 1 + k; < p < H — ). Next lemma yields 
Sobolev estimates for f n in terms of the Sobolev norm ll/xllvK^fC)- 

Lemma 5.2. For any 1 < p < 2 there exists a constant C = C(p) such that 

ll/nlk-(c)<c(p)r K n - 1 (M(p)) n J 

for any n = 1, 2, .... 

Proof. To prove the Lemma, we will use Remark l3.2l for the complex dilata- 
tion (j, and an arbitrary function g € W a ' p . Then it holds that 

\\D a (M)\\ LP{C ) < K\\D a g\\ L v m + CT \\g\\ *, (5.9) 

L'-'-p(C) 

for some positive constant C = C(p) > 1. First of all, we study the L p norm 
of f n . Recalling that \x is compactly supported inside of B, we first see that 

ll/n||iP(C) = ||/n||i>(D) < K \\T n fn-l\\LP(B)- 

Next, ([53]) yields that, 

\\D a fn\\LP(C) = \\D a {^Tnfn-l)\\LP{C) 

< CTq T n / n _l|| 2p_ -f K \\D a T n f n _i\\ L p(n\ 

L?=p(C) 

Hence, for any n > 1, 

||/n||w*.P(C) = ll/n||LP(C) + \\D a f n \\LP(C) 

<CTo|T n / n _i 2p_ + K\\T n f n -i\\ Wa , P t c \ 

To control the first term above, we see that 

H^n/n-lH 2p_ <-D||/n-l|| Jv_ < (Dk) ||T n _ 1 / n _ 2 || 2p_ 



< (Dk)- 1 ||Ti/o|| 2p <(£>k) t 
and for the second , if n > 1 

||^n/n-l||vy a -P(C) — -^ ll/«-l|lvK a -P(C)- 
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lp 2 



If we denote X n = \\fn\\w a >p(C) then we have just seen that 

X n <C 1 (KD) n + (KB)X n _ 1 (5.10) 

whenever n > 1, and where C\ = CTq\B>\p 2 . For n = 1 we proceed dif- 
ferently. Since both Ti/o and D a Tifo belong to L 2 (C), we can use Holder's 
inequality to get 

-— - 
1 1 ^1 /oil LP (B) + II -^"^l /o|| LP (D) ^ (II^i/o||l 2 (C) + II-O^i/oHl^C)) l^l p 2 



i i 

2 



< IBIp aS||/nllwa,2m = IBIp * BTo 



\Tifo\\w a ' 2 (C) l^l p 

fo\\w a - 2 . 

Thus 

Xi < C x kD + B \J])\p~^ T 

Thus, after recursively using (|5.10p . we end up with 

n— 1 -p 

X n < Ci K n Y, B j D n ^ + ( K B) n — |D|i~5 < ^ K «-i (^^ + D )n 

3=0 

where Ci = max{Ci,r |D|p~ 2 }. N ote finally that 

Ci<C<j>)To, 

which yields the claim. D 

In particular, every function /„ of the Neumann series is compactly sup- 
ported and belongs to L P (C) for any p G (l,oo), and also to W a ' p (C) for 
any p < 2. 

Lemma 5.3. If h belongs to W a ' p (C) for some 1 < p < 2, then 

|/<OI<<M(oy<C( rt ^=H 

\i\>R J it- 

Proof. We wil use the characterization in terms of Bessel potentials of 
W a ' p (C). Since the Fourier transform maps continuously L P (C) into L 9 (C), 
we get that 

c a+\t\ 2 )%\h(o\) q dA(oy<c( P )\\h\\ a>p 
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Thus, a simple computation yields 



\h(0\ 9 dA(O) < ( [ 



(i + l^l 2 ) 2 



a \ q 



\HO\ q dA(0 



< 



1 






<C(p) 



a,p 



R a 



and the result follows. 

Proof of Lemma I5.il We use the Neumann series 



*1. \ • 



ra+1 



-(n+l)fc Jn 



D 



(5.11) 



introduced before. Then, take g = Xm=o (f Ae_fc //T) (lAe-fc^) and 
/i = 5jV ~~ 5- I n this way, properties 1 and 2 follow easily from the general 
theory of the Beltrami equation, since 



Xe-kl^T 



" 'k 



\e_kn 

< k\\T\\ 

< {n\\T\\ s ) n \\v\\ s = {K\\T\\ s ) n k 



-\e_ k nT\ i-Xe-k/j, 



For the proof of 3, we must use the regularity of \i. Use l5.1l1 to write g(z,k) 

n+l 



YZ=o G n{k,z) where G n {z,k) = (|Aj e„ (n+1)fc / n . Then, Lemma 

can be applied to f n . The Fourier transform of G n (z, k) (with respect to the 
z variable) reads as 



G n (Z,k) 



k \ n+1 

^Aj / n (f-( n + i)fc) 
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1 1 

1 "0 / [■ \ q 



Hence, for \k\ > R, we can use lemma IST31 to get 

( / \gfo k)\" dA(0) " < E ( / \G n (£, fc)|« dA(0 
\ J \Z\<R J ^\J\e\<R 

no / » 

= E / |/»tt-(n+l)A0|«dA(£) 
„=o V/KK* 

no / „ \ j 

=E / i^or^o 

n=0 V^IC+(n+i)fe|<R / 

no / » \ i 

<E / |/n(C)| 9 ^(C) 

„=o \^ICI>(n+i)|fe|-fl / 

"0 

< cw e 



n ° II f II 

\\Jn\\a,p 



((n + l)\k\-R) a 

n=0 7 

where C(p) is the constant from Lemma 15, 31 Now, using Lemma E 

~ ^- no ( K M(p)) n 



J ■ m ,w*MQY<c M ±z ;^ 



H)< 



n=0 

r n ° i 

< c ^p)(«jf M) --ag ^- npF -^ ; 



and if we take \k\ > 2i£, then we finally get 

\At\<R J K \ k \ ^o( n + l 

<C(a,K,p)M(p) no 



y 

r i 



K \k\ a 

and the result follows. □ 

The Cauchy kernel is not in I? but it belongs locally to W e,p for 1 < p < 2, 
e < — j£_ Thus we can work with a mollification of it which is perfectly 
controlled. However we need to choose carefully the mollification kernel (see 

E2I vol l &v.i). 



Lemma 5.4. There exists a C* > such that for any N > 0, there exists a 
C°° function <pN in C having the following properties: 

• < (J)n < 1; 0jv = 1 on ID and 4>n = on 2B. 
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• f(f> N = 1. 

• \D a (p N \ < (C*iV)l Q l /or any a£l? + with \a\ < N. 
Lemma 5.5. Let $(z,w) = j^ and 1 < p < 2. 

(a) ||$(-,z)|| LP(D) < C(p) for all z G C. 

(b) $(•, z) € W £ ' p for e < —^ uniformly in z. 

(c) For any N > 0, there exists a mollification ^s,N such that 

\\H;z) ~ ® 5,n (-,2)11 lp (B) <C(e,p)5 e 
whenever z£C and e < —^-- 

(d) ||$5,jv||l2(c) blows up as a power of 5, i.e. 

\\^8,n(-,z)\\ L 2 {c) <C{p)5 l ~v 

(e) For each R > $ and m > 0, there exists a universal constant C* and 
C = C(p) such that for any m < N 

\\^N(;z)\\ L2m > R) < C(p)(C*N) m 5 l -l (5R)- m 

Proof. Claims (a) and (b) follow by the compactness of the support and 
Lemma 13.11 Now define 

Claim (c) follows from the fact that since p < 2, W e,p C B P ' ()3.2j) . Namely, 
||*z(0 -^8,n{z,-)\\lp < / Up{9 z )(w)<l>s(w)dw 



L2(C) 



< W*zU, p , 2 J(Mw)) 2 W 2+£2 )* <8 e (J<t>\y)\y\ 2+e2 f* <<f 



For claim (d), using Plancherel, Holder, Hausdorff- Young inequalities and 
(a), we obtain, for 1/p — 1/q = 1/2, that 

WHnWv < ||$*IMII^O)IU« < C5 l -'p. 

For the last claim, write again 

\\®6,N\\L 2 (\t\>Ro) < \\$z\\lp\\\(I>n(8Z)\\l<i(\Z\>Ro) 

< \\**\\lfi6 1 - 2/p \\MZ)\\L<(.\t\>SBo) 
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Now 



, (a + H2T 



\\<t>N{£)\\Li{\t\>5Ro) < III |7T^ \4>n{Q\\li{\Z\>5Ro) 

<(5Ro)- m \\ e ^i5^(e)iiiM<(^or m E ^ii^^Hl,' 

~— ' a! ^^ a! 

|a|=m |a|=m 

< (^ ) _m Y ^{C*N) m < (5R )~ m (2C*N) m 
* — ' a! 

|a|=m 



for m < N from where (d) follows. □ 

Now we combine the above estimates to obtain the precise decay for the 
solutions to the linear equation. 

Proposition 5.6. Assume that u G W a ' 2 (C) is a Beltrami coefficient, with 
compact support inside of D, such that ||/u||oo < K an d ||A*||w a < 2 (C) — ^o- 
For each A G 9B and eac/i fcdC, let tp = ip\(z,k) be the quasiconformal 
mapping satisfying 

— k 

dip\(z,k) = -Xe- k (z)^(z)dipx(z,k) (5.12) 

and normalized by 

ip\(z, k) — z = 0(l/z), z — > 00. 

There exists positive constants C = C(k) and b = b(n) such that 

\Mz,k)-z\<j^ 

for every z, k G C and every A G <9D. 

Proof. Let b > a constant to be defined , and let no G N. As in [T5], we 
can represent 

Mz,k)-z = C f^^-dA(w) 
Jo w-z 

= C ®(w,z)(g(w,k) +h(w,k))dA(w) 

Jc 

with g = g\(z, k) and h = h\{z, k) as in Lemma I5TT1 

Recall that we have control on g for low frequences by property 3 in 
Lemma l5.lt whereas h will be controlled by the ellipticity. It is also conve- 
nient to consider the mollification ^s,N of $ given in lemma 1531 for N to be 
chosen along the proof. We will therefore estimate the following four terms 
separately. The first three are dealt with by the usual ellipticity theory and 
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the Sobolev regularity of the Cauchy Kernel. Hence the estimates will de- 
pend on a suitable exponent s = s(k). It is in the last term where a,p will 
appear. Then we will chose the exponent p = 4/3 which will yield better 
constants. 



1= / ®(w,z)h(w)dA(w), 

Jo 

11= / ($(>, z) - <5>s,n(w, z))g dA(w), 

Jo 

111= / ^(Z,z)g(£,k)dA(0, 



IV= / $ StN (t;,z)g({i,k)dA(0 
J\Z\>R 

I:The tail Fix s = s(k) such that k IITIL < 1. Then we have 



$(w, z) h(w) dA(w) 



< \\$(;Z) 



<C(k,s) (k\\T\\ s ) 



<<() 



log(|fe|) 



-log(«l|r||.) 

C(K, S )(K\\T\\ s ) no <\k\- b 

in < lifer 6 . 



C(«)(l + 61og(|fc|) 



since by Lemma [531 (a), the norm ||$(-,z)|| *> does not depend on z. 
Take now, 

no > C(k, s) + b 

so that, 

and hence 



(5.13) 
(5.14) 



II: The error of mollification. We will use Lemma 15.51 with exponent 
1 < s' < 2 with 1 + jr = 1. Thus 0<e<l-fand0<5< \k\^r . Then it 
follows from Lemma 15.51 (c) and Lemma ET] that 

|H| < \\g\\ L s m \m;z) - ^ 5>N (;z)\\ L ^ x < C(K,s,e)5^ < \k\- b . 
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Ill: The mollification at high frequencies . We now choose the 
optimal value of N. We use Plancherel's Theorem and Lemma 15.51 (e) with 
m = N (assuming RS > 1), to get 



W£,*)?(£,&)cM(£) 



t\>R 



< lbllL2(C) \\^8,n(-,z)\\l^(\^>R) 

< C{s,k)(C*N) n 5^- 1 (SR)- n < \k\- 



(5.15) 

Let us plug in the value of 5 and choose N to obtain the optimal value of 
R. Namely first 5~~ ~ |A:| 2b . Thus we obtain that 



N\i,\3b+ — 



or 



iT > (CA^) iV |A:| 



R > C\k\^N\k\^ . 



With the optimal N = [3b\og(k)] + 1] we get the condition 

U>C|Jfe|Slog(|ife|). 
Imposing b < e we obtain that for large \k\ it is enough to take 



(5.16) 



IV: The mollification at low frequencies. The final term is the crucial 
one. Take 1 < p < 2, and q = — ^y. Then 



&<R 



g(Z,k)$6,N(Z,z)dA(w) 



< 



\t\<R 



m,k)\*dA(£)\ \\^ N (; Z )\\ LPm<R) 



For | A; I > 2R we can use Lemma 15. II and obtain 

m,k)\UA(0) <C(a, K ,p)M(p) no 
\£\<R J 



k\ a 



At the same time, the other factor is bounded with the help of Lemma 15.5 
(d), which is allowed since p < 2. More precisely, we have 



\$S,n(-,z)\\lp(\z\<r) 



\<R I 



< C{p) Rp 



±-i 



\^s,N^,z)\ 2 dA(0 



. '\t\<R 

<C(p)rI- 1 \\^ 5:N (;z)\\ L 2 {c) 

R\?~ 1 . „,(f )(!-i) 



<C(p) ( - ) < \k\ y ~ !{ p 
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Here we have inserted the values of R and 5 from II and III. Thus, whenever 
\k\ > 2R > y we have 



g(£,k)$ SiN (£,z)dA(w) 

\e\<R 



<C(a,K,p)^M(py- ' '~ ul 



9_ J\/Tf^\ n IZ,|(t)(p _1 ) 



Now since ||T||x,p < C(p — 1) it follows that the best choice is p = 4/3. 
Inserting this and the value of no from (|5.13|) in the previous equation, 

C(a, k, P ) T -!- \k\ c W b \k\i< C(a, K,p) r l^l 60 ^)^ 1 — 

Finally we want that (IV) is controlled by k as well. Since e = e{n) < 1 
and we already asked b < e, we end up getting that it suffices that 



f ea i ea 
6 < nun (-,*) = -. 



C" J C 
Here C = C(k) > 1 and we have use that a < 1. The proof is concluded. □ 

5.2 Estimates for the nonlinear equation 

Now that the behavior at k — > oo of the solutions to the linearized equation 
(|5.12p is known, it is time to study the behavior of the complex geometric 
optics solutions. 

Theorem 5.7. Let pi € W a ' 2 (C) be a Beltrami coefficient, real valued, com- 
pactly supported in jD, such that \\fi\\oo < ~kTT an< ^ Wl^Ww^^OC) — ^o- Let 
cj) = 4>\(z, k) be the solution to 

_ ~fc 

d(j>x(z,k) = --\fx(z)e-k(,<j)\(z,k))d(f>\(z,k) 

<ft\(z,k) — z = 0(l/\z\) as \z\ —> oo. 
There exists constants C = C(K) > and b = b{K) such that 

nv K 

\4>x{z,k)-z\< ° 



\k\ ha 

for every z € C, k E C and A £ dD. 

Proof. Since the estimate we look for is uniform in z and A, it suffices to 
show equivalent decay for the inverse mapping ip\ = 4>^ . But ip\ is the 
only quasiconformal mapping on the plane that satisfies both the equation 

— k 

dipx(z,k) = T \e-k(z)/j,(iJ)\(z, k))dip\(z,k) 
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(compare with (J5.12J1 ) and the condition i()\(z, k) — z = 0(l/\z\) as \z\ — ► oo. 
Then, we just need to show that the coefficient 

satisfies the assumptions of Proposition 15.61 First, it is obvious that 

K -1 

and it is also obvious that supp(/U o ifr\(-, k)) CB (this follows from Koebe's 
2 Theorem). Then, it remains to prove that \x o ip^ £ W@ ,2 (C) for some (3 € 
(0, 1). But this follows from Proposition 14.21 Indeed, since /i € W a,2 (C) n 
L°°(C), we have /i o Va € W^' 2 (C) with 

j_ 
I|^°^a|Iw^ 2 (C) ^ C r (f 

for any < (3 < j^, where C = C(a,/3,K). Note also that (5 behaves 
linearly as a function of a, with constant depending only on K. So the 
result follows. □ 

Remark 5.8. In the above result, the assumption supp(/u) € jB is not 
restrictive. Indeed, if supp(//) C D(0,R) for some R > then the function 
I^r(z) = /i(4i?z) defines a new Beltrami coefficient, compactly supported in 
|B, does not change the ellipticity bound, and 

One can then apply the previous Theorem to this coefficient /x/j and ob- 
tain estimates for the complex geometric optics solutions. But f/j, B (z,k) = 
ffj,(4Rz, -£g) and in fact if we represent these solutions as /^(z, k) = 
exp(ik(f>n(z,k)), then 



k+'fa'Ik 



W z > k ) = 777 <?V I 4jRz > 775 



so the estimates for 4>^ R coming from the previous theorem give similar 
estimates for </>„, modulo a power of i?. 

Now as discovered in |11] the unimodular complex parameter A allows to 
push the decay estimates to complex geometric optics solutions to the 7- 
harmonic equation. As always, given a real Beltrami coefficient v we denote 
by fi/(z,k) = e lkz M„(z,k) the complex geometric optics solutions to df = 
udf. 

Theorem 5.9. Let \i be as in Theorem \5.T\ and define 

u = Re(/ M ) + i Im(/_ /U ). 

There exist a function e = e(z, k) and positive constants C = C(K) and 
b = b(K) such that 
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(a) u(z,k) =e ifc (* +e (*' fc ». 



(b) \e(z,k)\< 



CT K 

- ! - for each z, k E C. 



1*1 



b a 



Further, a similar estimate holds for u = Re(/_ /x ) + i Im(/«) 
Proof. A calculation shows that it may be rewritten as 



u = fu 






1 + 



7/i ■/- 

/„ + /- 



Thus, the Theorem will follow if we find a function e(z, k) such that 



[e(z,fc)| < ^r and 



< 1-e 



fce(z,fe)| 



Following [TTJ Lemma 8.2], it suffices to see that 



inf 
t 



J» J-H i e «t 



/«+/- 



>e 



fce(z,fc)| 



For this, define $t(z, k) = e 2 (/^ cos(i/2) + i/-^ sin(£/2)). It follows eas- 
ily that for each fixed fc, 

J|e-^*$t(*,fc)-l|=0(l/*) as|z|^oo 

Thus, by uniqueness in Theorem 12 .2\ $t is nothing but the complex geomet- 
ric optics solution <J> t = f\^ with A = e . But then 



fu + f- 



V 



JfJ, "T /— (U //i 1 + 



A7_ 



M„ 



On the other hand, from Theorem 15.71 we get that 

e -|fc «(*,fc)| < |M M (z,A;)| = | e ifc ^"(*'*)-*)| < e |fce(2 ' fc)l 



where |e(z, fc)| < -77T&- and 



|fe|' 



-2\ke(z,k)\ < 



l/u(*,fc)l 



< e 2|fce(*,fc)| 



uniformly for A € e?D. Finally, by Theorem 12.21 we a l so have Re 
0, so that the result follows. 



M„ 



> 



D 
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fm ||wi,(i+c«e)2f D i < c (1 + r ) K 



6 Proof of Theorem 11.11 

In order to get stability from A 7 to /i, we will need the stability result for the 
complex geometric optics solutions in terms of a Sobolev norm. It comes as 
an interpolating consequence of the L°° stability result given at Theorem l5.9l 
and of the regularity of the solutions to a Beltrami equation with Sobolev 
coefficients (see Theorem I4.6p , 

Theorem 6.1. Let H\,\ii be Beltrami coefficients, compactly supported in 
D, such that \\idj\\ < J^fj and \\fij\\wa : 2(Q < Tq. Let f^. denote the complex 
geometric optics solutions to df flj = fijdf^. Then, for each 6 € (0, -^) we 
have 

1 -ba 2 

log- 
P 

for come constants C = C(\k\, K) > and b = b(K) > 0. In particular, the 

same bound holds with the W 1,2 (J$)-norm. 

Proof. The subexponential growth obtained in Theorem 15.91 entitled us to 
apply Theorem 12.41 B to the solutions u 7i . Since they are equivalent to the 
corresponding /„ we achieve the estimate 

i 

\\Ud-,k)-Ud-M\L^ 0) <j^^ 

for some positive constants C = C(k, K) and b = b(K). On the other hand, 
from Theorem 14.61 f° r every 9 G (0, -^) we have 

II/a*(->*0 ~ f^(^ k )\\w^+-<>,2(B) = \\D 1+9a (/ M (-,fc) - / M2 (-,fc)) |[ L 2 (D) 

<Ce c ^ (l + \k\) (T + \k\ a ) e . 

As in Theorem 14.61 here C may depend on K. Let ip € C°°(C) be a cut-off 
function, compactly supported in 2D), such that <p\n = X\o- Then the above 
estimates imply that 



< 



CY 



\Kfm(;k)-f„ 2 (;k))<p\\L*c iC ) - |l og ( p )|6a 

||(/ w (•, k) - U 2 (; k)) <p\\ Wl+a e, 2{c) < C e c ^ (1 + \k\) (To + \k\ a ) e 

where, as usually, || • ||^ s , P (<n denotes the homogeneous Sobolev norm. Now, 
an interpolation argument shows that for each < (3 < 1 we have 

\\<P(UA-,k)-U 2 (-,k))\\ w(1+ae)Pt 2^ 

1-/3 

< C e c ^ (1 + \k\f (T + \k\ a f e ^f 



l 0g (p)|6a(l-/?) 
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where C = C(K). In particular, for (3 = j^q we get that 
Wffii-f^Ww't^+'x^iO) — II (/mi ~ f^) i P\\w 1 '( 1 + ae ') 2 (C) 

< C e°\ k \ (1 + |fc|) (1 + r )^ +T r^ (9 ~^ ) - 



ba 2 6 

log(p)|i+«fl 



Here the sharp modulus of continuity is obtained when the logarithm has 
the bigger exponent, which is given for 9 = j^. Thus, we end up with 

II /w -/wllw-i-a+a^an ^(W) (i + r )^ — rr 

|log(p)|* a 

as claimed. D 

It just remains to see how the previous estimate drives us to the final stability 
bounds for the Beltrami coefficients (and therefore for the conductivities). 
To do this, the following interpolation Lemma will be needed. Note that it 
includes L p spaces with p < 1. 

Lemma 6.2 (Interpolation). Let < po < 2 and 2 < p\ < oo. Let 9 be 

such that 

1 _ _#_ 1-0 

2 po p\ 
Then 

II/IIl 2 <II/I|£,oII/II^? 

for any f G D>° D L Pl . 



Proof. The proof is adapted for the usual Riesz method for interpolation 
with a little extra care when po < 1. We choose r < po and define exponents 
QOiQii Q2 such that 

111 111 111 

- = — + - - = - + — _ = _ + _ 

r po q r 2 q 2 r p\ q x 

For z = x + iy in the strip < y < 1, we define the analytic function 

G(z) = \ g \^ + ^f r 

\g\ 

Notice that \G(iy)\ qi = \g\ q2 , |G(1 +iy)|* J = \g\ q \ and |G(0 + ij/)| = \g\. 
Now we introduce the function 



!{*)=[ \f\ r \G{z)\ 
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We can estimate its values at the boundary of the strip, 



\im<\\fhpi\\G(iy)\\w= H/ll 



LPi 



(/ 



\<n 



\I(l + iy)\<\\fhpo\\G(iy)\\ m =\\f\\LPo 



\9 



'12 



Then we apply Phragmen-Lindelof theorem to the function I(z) obtaining 
that 



I(e + iy)< [\\f\\ LP1 



92 



91 \ 



l/ll 



LPO 



'12 



1-6 



But 1(6 + iy) = ||/<7||v, so the result follows. 



D 



We are finally led to obtain the desired stability in L 2 norm of the Beltrami 
coefficients. 

Corollary 6.3 (Proof of Theorem II .ip . Let ii\,[i 2 be Beltrami coefficients, 
compactly supported in B, such that \\fij\\ < ^Ej an d llMjllw a < 2 (C) ^ ^V 
There exists constants b = b(K) > and C = C(a, K) > such that 

-be? 



Imi -M2||l2(d) < c (i + r )^ 



log' 



where p = II Ai — Aoll i 



i 



Proof. Denote by /j the complex geometric optics solution A, 4 of 9/ = /ij 5/ 
with fe = 1. Then, 



I Ml — A*2| 



dhdh-df 2 dh 

dfi df 2 


= 


-oh (a/i - a/ 2 ) + (a/i - 0/2)0/1 


dfi df 2 



, \dfi-df 2 \ , 



\dfi-df 2 \ ^\Dh-Df 2 \ 



< 2 



\df 2 \ |9/ 2 | |5/ 2 | 

because \Dfj\ = \dfj\ + |0/j|. Therefore, for any s > 

£>/i - £>/ 2 



[Ml - M2||l s (D) < 2 



9/2 



L s ( 



Now, let p e (0, jfc) and G (0, l/K). Then put A - ^j+jg -r - 
application of Holder's inequality gives us that 

Dfi - Df 2 



+ h An 



df 2 



<C(a,6)\\Df l -Df 2 \\ L2{1+a e HD) 



< C(a,0) ||/i - /2 [1^-1,3(1+019) ( 



1 



df 2 
1 



Lp( 



df 2 



LP( 
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Now, using Lemma 14.71 and Theorem 16,11 we obtain the estimate 

1 



i 



I Mi - frh-cn) <c(i + t )k 



log 

p 



-bo? 



where C > depends on a, 6, and K, and b > depends on i'f . Finally, 
if s > 2 then we are done, since /Uj are compactly supported in D. But in 
general we only know < ^ < s so one could well have s < 1. In this 
case, in order to get L? estimates only interpolation between L s and L°° is 
needed, as in Lemma 16.21 Namely, 



2-i 

2 II,, ,, II 2 



llMl _ M2||l 2 (B) < II Ml ~ M2 1 1 £, S (][J) ||Mi _ M2|| L . 

and now the stability estimate looks like 

- ba 



log- 

P 



IImi -M2||l2(d) < c (i + r )^ 

where the constants may have changed. □ 



References 

[1] D.R. Adams, Sobolev spaces, Academic Press, New York, 1975. 

[2] D.R. Adams, L.I. Hedberg, Function spaces and potential theory, 
Springer (1996). 

[3] L. V. Ahlfors, Lectures on quasiconformal mappings. Manuscript 
prepared with the assistance of Clifford J. Earle, Jr. Van Nostrand 
Mathematical Studies, No. 10 D. Van Nostrand Co., Inc., Toronto, 
Ont.-New York-London 1966 

[4] G. Alessandrini, Stable determination of conductivity by boundary 
measurements. Appl. Anal. 27 (1988), no. 1-3, 153-172. 

[5] G. Alessandrini, Singular solutions of elliptic equations and the de- 
termination of conductivity by boundary measuraments, J. Diff. Eq. 84, 
(1990), 252-272. 

[6] G. Alessandrini, S.Vesella Lipschitz stability for the inverse con- 
ductivity problem. Adv. in Appl. Math. 35 (2005), no. 2, 207-241. 

[7] K. Astala, Area distortion of quasiconformal mappings, Acta Math. 
173 (1994), 37-60. 

[8] K. Astala, D. Faraco, L. Szekelyhidi Jr. Convex integration and 
the LP theory of elliptic equations. Preprint (2004). 

47 



[9] K. Astala, T. Iwaniec, G. Martin, Quasiconformal mappings and 
PDE in the plane. Monograph in preparation. 

[10] K. Astala, T. Iwaniec, E. Saksman, Beltrami operators in the 
plane, Duke Math. J. 107 (2001), no.l, 27-56. 

[11] K. Astala, L. Paivarinta, Calderon Inverse Conductivity problem 
in plane. Ann. of Math. (2) 163 (2006), no. 1, 265-299. 

[12] K. Astala, L. Paivarinta, A boundary integral equation for 
Calderon Inverse conductivity problem. Collect. Math. 2006, Vol. Ex- 
tra, 127-139. . 

[13] K. Astala, M. Lassas, L. Paivarinta, Calderon's inverse problem 
for anisotropic conductivity in the plane. Comm. Partial Differential 
Equations 30 (2005), no. 1-3, 207-224. 

[14] J. A. Barcelo, T. Barcelo, A. Ruiz, Stability ot the inverse con- 
ductivity problem in the plane for less regular conductivities, J. of Diff. 
Eq. 173, 2001, 231-270. 

[15] T. Barcelo, D. Faraco, A. Ruiz, Stability of Calderon inverse 
conductivity problem in the plane, J. Math. Pur. Appl., 2007; 88 (6), 
522-556. 

[16] R. Beals, R. Coifman, Multidimensional inverse scattering and non 
linear partial differential equations. Pseudodifferential Operators and 
Applications (F. Treves ed.) Proc Sympos. pure math. Vol 43, pp 45- 
70. AMS Providence 1985. 

[17] J. Bergh, J. Lofstrom, Interpolation spaces. An introduction, 
Springer New York 1976. 

[18] R. Brown, Global uniqueness in the impedance-imaging problem for 
less regular conductivities, SIAM J. Math. Anal. 27 (1996), no. 4, 1049- 
1056. 

[19] R. Brown, Recovering the conductivity at the boundary from the 
Dirichlet to Neumann map: A pointwise result, J. Inverse Ill-Posed 
Probl. 9 (2001), no. 6, 567-574. 

[20] R. Brown, Estimates for the scattering map associated with a two- 
dimensional first- order system, J. Nonlinear Sci. 11 (2001), no. 6, 459- 
471. 

[21] R. Brown, R. Torres, Uniqueness in the inverse conductivity prob- 
lem for conductivities with 3/2 derivatives in L p . J. Fourier Analysis 
and Appl., 9, 6 (2003) 563-574. 



18 



[22] R. BROWN, G. UHLMANN, Uniqueness in the inverse conductivity 
problem, for nonsmooth conductivities in two dimensions, Comm. in 
PDE, 22, (1997), 1009-1027. 

[23] A. P. Calderon, On an inverse boundary value problem, Seminar on 
Numerical Analysis and its Applications to Continuum Physics, Soc. 
Brasileira de Matematica, Rio de Janeiro, (1980), 65-73. 

[24] A. Clop, D. Faraco, J. Mateu, J. Orobitg, X. Zhong, Beltrami 
equations with coefficient in the Sobolev space W l,p , Publ. Mat. to 
appear. 

[25] A. Clop, X. Tolsa, Analytic capacity and quasiconformal mappings 
with W 1 ' 2 Beltrami coefficient, Math. Res. Lett, to appear. 

[26] J. Duoandikoetxea, Fourier analysis, Graduate Studies in Mathe- 
matics 29, American Mathematical Society, Providence (2000). 

[27] D. FARACO, Milton's conjecture on the regularity of solutions to 
isotropic equations, Ann. Inst. H. Poincare Anal. Non Lineaire 20 
(2003), no. 5, 889-909. 

[28] F. Gehring. Characteristic properties of quasidisks. Seminaire de 
Mathematiques Superieures [Seminar on Higher Mathematics], 84. 
Presses de l'Universite de Montreal, Montreal, Quebec, 1982. 97 

[29] E. Hernandez, G. Weiss, A first course on wavelets. With a fore- 
word by Yves Meyer. Studies in Advanced Mathematics. CRC Press, 
Boca Raton, FL, 1996. 

[30] T. Iwaniec, G.Martin, Geometric function theory and non-linear 
analysis, Oxford Mathematical Monographs, The Clarendon Press Ox- 
ford University Press, New York, 2001. 

[31] C.E. Kenig, A. Ponce, L. Vega On the generalized Benjamin-Ono 
equation. Trans. Amer. Math. Soc. 342 (1994), no. 1, 155-172 

[32] R. V. Kohn, M. Vogelius, Identification of an unknown conduc- 
tivity by means of measurements at the boundary, Inverse problems 
(New York, 1983), 113-123, SIAM-AMS Proa, 14, Amer. Math. Soc, 
Providence, RI, 1984. 

[33] R. V. Kohn, M. Vogelius, Determining conductivity by boundary 
measurements. II. Interior results. Comm. Pure Appl. Math. 38 (1985), 
no. 5, 643-667. 

[34] K. Knudsen, On the inverse conductivity problem, Ph. D. Thesis. 
Aalborg University 2002. 



49 



[35] M.Lassas, J.L.Mueller, S.Siltanen, Mapping properties of the 
nonlinear Fourier transform in dimension two, To appear in Comm. 
Partial Differential Equations. 

[36] L. Liu, Stability estimates for the two-dimensional inverse conductivity 
problem, Ph. D. Thesis, Dep. of Mathematics, University of Rochester, 
New York, (1997). 

[37] O. Lehto, K. I.Virtanen, Quasiconformal mappings in the plane, 
second ed. Springer- Verlag, New York, 1973. Translated from the Ger- 
man by K. W. Lucas, Die Grundlehren der mathematischen Wis- 
senschaften, Band 126. 

[38] N. Mandache, Exponential instability in an inverse problem for the 
Schrodinger equation, Inverse Problems 17, 2001. 1435-1444. 

[39] A. Nachman, Global uniqueness for a two dimensional inverse bound- 
ary problem, Ann. of Math. 143 (1995) 71-96. 

[40] C. Perez, Banach function spaces and the two-weight problem for 
maximal functions, Function spaces, differential operators and nonlin- 
ear analysis (Paseky nad Jizerou, 1995), 141-158, Prometheus, Prague, 
1996. 

[41] L. Paivarinta, A. Panchenko G. Uhlmann, Complex geometrical 
optics solutions for Lipschitz conductivities, Rev. Mat. Iberoamericana 
19 (2003), no. 1, 57-72. 

[42] H. Renelt, Elliptic systems and quasiconformal mappings. John Wi- 
ley and Sons. New York 1988. 

[43] H. M. Reimann, T. Rychener, Funktionen beschrdnkter mittlerer 
oszillation, Lect. Notes in Math. 487 Springer- Verlag, Berlin (1975). 

[44] M. Schechter, Principles of functional analysis. Student edition. 
Academic Press [Harcourt Brace Jovanovich, Publishers], New York- 
London, 1973. 

[45] E.M. Stein, Singular integrals and differentiability properties of func- 
tions, Princeton Math. Series 30. Princeton University Press. N.J. 
1970. 

[46] D. Swanson, Pointwise inequalities and approximation in fractional 
Sobolev spaces, Studia Math. 149 (2) (2002), 147-174. 

[47] J. Sylvester, G. Uhlmann, A global uniqueness theorem for an 
inverse boundary value problem, Annals of Math., 125, (1987), 153- 
169. 



50 



[48] J. Sylvester, G. Uhlmann, A uniqueness theorem for an inverse 
boundary value problem in electrical prospection, Comm. Pure and 
Appl. Math 39,(1986) 91-112. 

[49] J. Sylvester, G. Uhlmann, Inverse boundary value problems at 
the boundary- Continuous dependence, Comm. Pure Appl. Math: 41 
(1988) 197- 221. 

[50] L.Tartar Lecture notes in homogenization and optimal design 
http: //www. math. emu. edu~nwOz/CIMECNA.dvi 

[51] H. Triebel, Theory of function spaces, Monographs in Mathematics, 
78. Birkhauser Verlag, Basel, 1983. 

[52] F. Treves, Introduction to Pseudo Differential and Fourier Integral 
Operators, University Series in Mathematics, Plenum Publ. Co. 1981. 

[53] I. N. Vekua, Generalized analytic functions, Pergamon Press, 
London-Paris- Frankfurt; Addison- Wesley Publishing Co., Inc., Read- 
ing, Mass. 1962. 

[54] W. P. Ziemer, Weakly Differentiate Functions: Sobolev Spaces and 
Functions of Bounded Variation, Springer Verlag. 



A. Clop 

www.mat.uab.cat/~albertcp 

Department of Mathematics and Statistics, 

P.O.Box 35 (MaD) 

FI-40014 University of Jyvaskyla 

Finland 

D. Faraco, A. Ruiz 
Departmento de Matematicas 
Universidad Autonoma de Madrid 
Campus de Cantoblanco, s/n 
28049-Madrid 
Spain 



51 



